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ABSTRACT 

Understanding the epoch of reionization and the properties of the first galaxies rep- 
resents an important goal for modern cosmology. The structure of reionization, and 
hence the observed power spectrum of redshifted 21cm fluctuations are known to be 
sensitive to the astrophysical properties of the galaxies that drove reionization. Thus, 
detailed measurements of the 21cm power spectrum and its evolution could lead to 
measurements of the properties of early galaxies that are otherwise inaccessible. In 
this paper, we make predictions for the ionised structure during reionization and the 
21cm power spectrum based on detailed models of galaxy formation. We combine the 
semi-analytic GALFORM model implemented within the Millennium-II dark matter 
simulation, with a semi-numerical scheme to describe the resulting ionization struc- 
ture. Semi-analytic models based on the Millennium-II simulation follow the properties 
of galaxies within halos of mass greater than ~ 1.4xlO 8 M at z > 6, corresponding 
to the faint sources thought to dominate reionization. Using these models we show 
that the details of SNe and radiative feedback affect the structure and distribution 
of ionised regions, and hence the slope and amplitude of the 21 cm power spectrum. 
These results indicate that forthcoming measurements of the 21cm power- spectrum 
could be used to uncover details of early galaxy formation. We find that the strength 
of SNe feedback is the dominant effect governing the evolution of structure during 
reionization. In particular we show SNe feedback to be more important than radiative 
feedback, the presence of which we find does not influence either the total stellar mass 
or overall ionising photon budget. Thus, if SNe feedback is effective at suppressing 
star formation in high redshift galaxies, we find that photoionization feedback does 
not lead to self-regulation of the reionization process as has been thought. 

Key words: Cosmology: theory; diffuse radiation; dark ages, reionization, first stars; 
Galaxies: high-redshift 



1 INTRODUCTION 

In anticipation of forthcoming 21 cm observations of the 
epoch of reionization, a great deal of theoretical attention 
has focused on the prospects of measuring the 21 cm power 
spectrum. To this end, significant progress has been made 
in modeling the effect of galaxies on the reionization of the 
IGM. In large modern simulations, the most common ap- 
proach is to begin with an N-body code to generate a dis- 
tribution of halos (e.g . |Ciardi et aL |2003 [Sokasian et al.| 
[20031 |llieveraLl|2^7i [Zahn et al.||2007| |Trac fe Cen||2007| 
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Shin et al.|2008| |Iliev et al.|2008| |Trac et al.|2008| ). A sim- 
ple prescription is then used to relate halo mass to ionizing 
luminosity. Following this step, radiative transfer methods 
(most commonly ray-tracing algorithms) are employed to 
model the generation of ionized structure on large scales. 
The radiative transfer is normally run with lower resolution 
than the N-body code for computational efficiency. 



These simulations describe the generic features of reion- 
ization (e.g.|Iliev et al.|2007||Zahn et al.|2007||McQuinn et al. 
2007||Shin et aI.|2008||Lee et al.|2008||Croft fe Altay|2008|), 



confirming expectations from analytic models (e.g. Furlan 



etto et al.|2004b|a| |Wyithe fe Morales|2007| |Barkana|2 009 ) 
that large-scale, over-dense regions near sources are gener- 
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ally reionized first, and that massive galaxies tend to be 
surrounded by clustered sources that increase the size of HII 
regions. In addition, the simulations describe the structure 
of the HII regions, showing that they are generally aspher- 
ical (even where the sources are assumed to emit isotropi- 
cally). The growth of HII regions during reionization may 
also be influenced by radiative feedback in the form of sup- 
pression of galaxy formation below the cosmological Jeans 
mass within a heated IGM (Dijkstra et al.|[2004[), although 



Furthermore Raicevic et al. (2011) and Benson et al. (2006) 



the importance of this effect remains controversial (Mesinger 
& Dijkstra 2008). Suppression of low mass galaxy-formation 



delays and extends the reionization process, which though 
started by low mass galaxies, must then be completed by 
relatively massive galaxies ( |Iliev et aL 2007). 

An important outcome from the large cosmological vol- 
umes attained by modern numerical simulations has been 
the prediction of 21 cm signals that will be observable us- 



ing forthcoming low frequency arrays (e.g. Mellema et al. 
2006( Lidz et al. 2008). The most generic features of 21 cm 



power spectrum modeling were elucidated by |Lidz et al.| 
(2008), who show examples of its evolution. On scales of 
k rsj O.lMpc -1 the amplitude and the slope of the power 
spectrum vary in a non-monotonic way relative to the ex- 
pected shape in the absence of ionization structure. Thus, 
measurement of the 21cm power spectrum will provide the 
first clues regarding the clustering of ionizing sources dur- 
ing reionization. In particular, Lidz et al. (12008) illustrate 



that the slope and amplitude of the 21 cm power spectrum 
vary considerably among different models at a given ion- 
ization fraction. However they also find that the behavior 
with ionization fraction across the different models is rel- 
atively generic. In particular, the amplitude of the 21 cm 
power spectrum reaches a maximum close to the epoch when 
~ 50% of the volume of the IGM is ionized, while its slope 
is found to flatten with increasing ionization fraction. |Lidz| 
|et al.| ( |2008| argue that first generation low frequency ra- 
dio telescopes like the Murchison Widefield Arrajjjand the 
Low Frequency ArrajFlwill have sufficient sensitivity to mea- 
sure the redshift evolution in the slope and amplitude of the 
21 cm power spectrum. 

One of the main limitations in modelling of reionization 
is the physics of the ionizing sources. Most studies have used 
very simple prescriptions to assign ionizing luminosities to 
dark matter halos. It has been shown that it is then possible 
to constrain the parameters for these simple prescriptions. 
However an important open question is the degree to which 
the important astrophysics governing formation and evolu- 
tion of high redshift galaxies is accessible via observations 
of the 21cm power spectrum. A few studies have previously 
addressed the issue of realistic modelling of high redshift 



galaxies. For example, Raicevic et al. (2011) (see also Ben 
son et alT]|2006 Lacey et al.||2011 ) used the semi- analytical 
galaxy formation code, GALF ORM (Cole et al.|2000||Baugir 



have studied the effect of SN feedback on the global ionizing 
photon budget and global ionization. However these studies 
were restricted to the global evolution, and do not address 
the ionization structure. 

In this paper we combine detailed models of high red- 
shift galaxy formation using GALFORM with calculations 
of the spatial dependence of reionization, and predict the re- 
sulting redshifted 21cm power spectrum. We begin in §|2]and 
§ H] by describing the implementation of GALFORM, and 
our method for modelling the ionization structure. Then, in 
§[4] we present ionisation maps for different galaxy formation 
models, including the effect of SNe and radiative feedback. 
We discuss the ionising photon budget as a function of halo 
circular velocity in § [5] and then present a discussion of the 
dependence of the 21cm power spectrum on the galaxy for- 
mation model in § [6] We finish with some conclusions in 
§0 



2 THE MODEL 

In this section we introduce the theoretical galaxy formation 
modelling used in our analysis. In § |2.1| we briefly review 
GALFORM. We then describe the implementations of SNe, 
AGN and photoionization feedback processes in § |2.2| 



2.1 The GALFORM galaxy formation model 

The formation and evolution of galaxy properties are com- 
puted within the ACDM structure formation framework us- 
ing the semi-analytical model GALFORM. GALFORM in- 
cludes a range of processes that are thought to be impor- 
tant for galaxy formation, including: (1) The gravitation- 
ally driven assembly of dark matter haloes; (2) The density 
and angular momentum profiles of dark matter and hot gas 
in haloes; (3) The radiative cooling of gas and its collapse 
to form centrifugally supported disks; (4) Star formation in 
disks; (5) Feedback processes, resulting from the injection of 
energy from Supernovae (SNe) and AGN heating; (6) Chem- 
ical enrichment of the interstellar medium (ISM) and hot 
halo gas which affect the gas cooling rate and the properties 
of the stellar populations in a galaxy; (7) The dynamical 
friction on orbits of satellite galaxies within a dark matter 
halo and their possible merger with the central galaxy; (8) 
The formation of galactic spheroids; (9) The spectrophoto- 
metric evolution of stellar populations; (10) The effect of 
dust extinction on galaxy luminosities and colours, and its 
dependence on the inclination of a galaxy; (11) The genera- 
tion of emission lines from interstellar gas ionized by young 
hot stars. 

A comprehensive overview of GALFORM can be found 
Cole et al.|(|2000|, with an updated discussion in the re- 



trees to evaluate the ionizing photon budget, finding that al- 
though galaxies should produce sufficient ionizing photons 
to complete reionization, most of the galaxies responsible 
would be below the detection threshold of current surveys. 



1 http://www.haystack.mit.edu/ast/arrays/mwa/ 

2 http://www.lofar.org/ 



et al.|2005| [Bower et al.|20 06), based on Monte-Carlo merger view by |Baugh| fl2006|). In this paper, we implement GAL- 



FORM within the Millennium-II cosmological N-body sim- 
ulation ( Boylan-Kolchin et al. 2009). However rather than 



using the halo merger trees presented in |Boylan-Kolchin| 
et al. (120091 , we base our study on the halo merger trees de- 



scribed in the study of Merson et al. (2012) which are better 



suited for the purposes of semi-analytic modelling. The sim- 
ulation has a cosmology including fractional mass and dark 
energy densities with values of Q m = 0.25 , Qh = 0.045 and 
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H a =0.75, a dimensionless Hubble constant of /i=0.73, and a 
power spectrum normalisation of as— 0.9. The particle mass 
of the simulation is 6.89xlO 6 h _1 M0 and we detect haloes 
down to 20 particles in the simulation box of side length 
L = 100/i _1 Mpc. 



2.2 Feedback processes 

Feedback processes during galaxy formation are very im- 
portant contributors to the shape of luminosity functions 



predicted by GALFORM (Cole et al. |2000; Benson et al. 



2002||Baugh et al.|2005||Bower et al.|2006||Kim et al.|2011 ) 
Three main feedback processes are implemented in GAL- 
FORM, which we discuss in turn below. 



2.2.1 SNe feedback 

SNe feedback on galaxy formation is implemented within 
GALFORM through the reheating and ejection of cold gas 
from galaxies via the equation 



AL 



eject 



M, 



(i) 



where tp is the instantaneous star formation rate. Here j3 is 
the efficiency of the feedback process, parameterized as 



13 = (Vkisk/K< 



(2) 



where Vhot has unit of kms -1 , ahot is a dimensionless ad- 
justable parameter which controls the strength of SN feed- 
back (see Cole et al. 2000| and Vdisk is the circular velocity of 
the galactic disk at the half-mass radius. SNe feedback sup- 
presses the formation of galaxies within small dark matter 
halos, and is required to reproduce the faint end of the ob- 
served galaxy luminosity function (e.g. |Norberg et aT]2002 



Blanton et al.|2001| ). [Bower et al.| ( |2006[ hereafter the Bow06 
model) adopted values of Vhot=485km/s and ahot =3.2. 



2.2.2 A GN feedback 

To reproduce the low number density of bright galaxies and 
the steep slope of the bright end of the galaxy luminos- 
ity function an additional feedback process is needed that 
operates in the high mass regime. For this reason, Bow06 
included AGN feedback which suppresses the cooling flows 
in massive haloes. The physical motivation for this lies in 
the fact that energy is known to be released from accre- 
tion of matter onto central supermassive black holes. Bow06 
modelled AGN feedback assuming halos to be in quasi- 
hydrostatic equilibrium in cases where the cooling time at 
the cooling radius, t CO oi(^cooi), exceeds a multiple of the free- 
fall time at the cooling radius, £ff(r CO oi), i-e. 



t'cool v cool J ^ 



C^cool 



-tff(r CO oi), 



(3) 



where a CO oi is an adjustable parameter whose value controls 
the strength of AGN feedback. The value of a CO oi in the 
Bow06 model is 0.58. An alternative approach was taken 
by Baugh et al. (2005, hereafter Bau05) who implemented 
superwind feedback from star formation with a top-heavy 
stellar initial mass function in order to understand the high 
luminosity end of the galaxy luminosity function. 



2.2.3 Photoionization feedback 

In the presence of a strong ionizing background, star for- 
mation in small galaxies is quenched owing to several 
physical processes, including the suppression of cooling by 
photo- heating ( Efstathiou||1992 ) , the higher IGM gas pres- 



sure (Gnedin 2000), and photo-heating (Hoeft et al. 



2006 



Okamoto et al. 2008). As a result the star formation rate 



density is suppressed within HII regions during reionization 
(see Crain et al.|2009 ), which may result in self-regulation of 
the reionization process ( Iliev et al.|2007 ). Based on Benson 
et al.| ( |2002[ ), GALFORM includes a prescription for sup- 
pressing the cooling of halo gas onto the galaxy when the 
IGM becomes globally ionized. In the standard implemen- 
tation this is assumed to occur at a particular redshift z cu t- 
It is assumed that the suppression of cooling occurs when 
the host halo's circular velocity lies below a threshold value, 
Kut, at redshift z < z C ut- Rather than a constant 2 cu t, in 
this paper we apply suppression when the cell in which the 
galaxy resides is fully ionized (see § 2.3). The adopted value 
is Kut = 30 km/s for each of the models (|Lacey et al.|2011| 



Lag os et al.|2012 ). We note that the value of Vcut = 30 km/s 
is used for the Bow06 model in this paper, rather than the 
original value Vcut = 50 km/s used in (Bo wer et al.J 2QQ6). 

2.2.4 Key differences of variant models 

In this paper, we use variants of the GALFORM galaxy 
formation model described in § |2.1| based on two main pub- 
lished implementations which we refer to as Bow06 (Bower 
et al.||2006| ) and Lagos ( |Lagos et aL|2012| ). The Bow06 and 
Lagos models assume a Kennicutt IMF, similar to that in the 
solar neighbourhood, for both quiescent star formation and 
starbursts. The Lagos model is similar to the Bow06 model 
but it has a new star formation law and has different pho- 
toionization feedback parameters. We study three variants 
of the Bow06 model, in addition to the published version. 
First, the NOSN model in which we remove SNe feedback 
from the Bow06 model by using a value of Vhot=0km/s. 
Second, the Bow06(no suppression) model in which we re- 
move photo- ionization by setting Vcut=0km/s, and third the 
NOSN (no suppression) model in which we remove both SNe 
and radiative feedback (i.e. Vhot = Vcut=0km/s). In Tablefl] 
we summarize the values of selected parameters for the dif- 
ferent models used in this paper. 



2.3 Modelling spatial dependence of radiative 
feedback in GALFORM 

Photoionization feedback from reionization is normally mod- 
elled in semi-analytic models (including GALFORM) using 
a single value of z cu t- However to investigate the effect of 
galaxy formation on ionisation structure during reioniza- 
tion we need to improve the photo-ionizing feedback using 
a spatially dependent value of z cu t that accounts for earlier 
suppression in regions of the IGM where HII regions first 
form. Broadly, the process is to run GALFORM to a par- 
ticular redshift snapshot, perform a calculation of ionisation 
structure, and then apply radiative suppression to subse- 
quent galaxy formation inside HII regions with halo circular 
velocities below Vcut- GALFORM then evolves the popula- 
tion of galaxies to the next snapshot, where the ionisation 
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Table 1. The values of selected parameters which are different in the models. The columns are as follows: (1) the name of the model, 
(2) the value of the photoionization parameter Vcut, (3) the SNe feedback parameter, Vhot, (4) the IMF of brown dwarfs T, and (5) 
comments giving model source or key differences from published models. 



V^t [kms 



Vhot [kms 



Comments 



)w06 


30 


485 


1 


LgOS 


30 


485 


1 


)w06(no suppression) 





485 


1 


3SN 


30 





4 


3SN(no suppression) 








4 



Bower et al. (2006), V cu t value change 

Lagos et al. (2012) 

Bower et al. (2006), No radiative suppression 

Bower et al. (2006), No SNe feedback 

Bower et al. (2006) 

No SNe feedback and No radiative suppression 



structure is recomputed. Details of the GALFORM specific 
implementation are provided in the Appendix. 



NOSN(no suppression) models. The value of T should be 
greater than unity by definition. 



2.4 The High redshift galaxy luminosity function 

It is important to consider how well the predicted galaxy 
population represents the galaxies observed to exist at high 
redshift. GALFORM models are calibrated to a wealth of 
data at low redshift. Previously, Lacey et al. (2011| used 



GALFORM with Monte Carlo merger trees (with no built- 
in mass resolution limit) to compare model predictions with 
observed properties of high-redshift galaxies at z ~ 3 — 10. 
Their modelling was successful in reproducing the high red- 
shift luminosity function, although with dependence on the 
model used. In our work we have utilised N-body merger 
trees extracted from the Millennium-II simulation in order 
to explicitly include the correlations between galaxy posi- 
tion and over-density in the IGM. Since the Millennium-II 
simulation resolves most of the galaxies responsible for reion- 
ization, the model predicts the correct star formation rate 
density. We show the UV luminosity functions at high red- 
shifts in Figure ^ for the Bow06 model and Lagos model 
based on the Millennium-II dark matter simulation merger 
trees. The UV luminosity function using the Millennium-II 
simulation merger trees is nearly identical to the Monte- 



Carlo merger trees for the Bow06 model (se e |Lacey et al. 
2011). The Bow06 model (solid lines in Fig. nl agrees well 
with the observational results, with the exception of an over- 
prediction of luminous galaxies. The Lagos model (dashed 
lines in Fig. [I]) shows better agreement with the data across 
all redshifts and luminosities. The Lagos model has a dif- 
ferent star formation law to the Bow06 model, and adopts 
a burst timescale that also gives better agreement with the 
UV luminosity function at z ~ 3 — 7 following the analysis in 
the Lacey et al. (2011[). Importantly, both models agree well 



for the faint galaxies thought to be responsible for reioniza- 
tion, indicating that our results should not be very sensitive 
to this choice. 

We also show the NOSN model. Simply removing the 
feedback strength of SNe (by setting Vhot = 0) results in a 
model that which greatly over predicts the number of galax- 
ies at all luminosities. In order to correct for this we there- 
fore modify the parameter in GALFORM which specifies 
the ratio between the sum of the mass in visible stars and 
brown dwarfs, and the mass in visible stars. This parameter 
(T) quantifies the assumption for the IMF of brown dwarfs 
(m < 0.1M©), which contribute mass but no light to stellar 
population. We adopt a value of T = 4 for the NOSN and 



SEMI-NUMERICAL SHEME TO 
CALCULATE THE EVOLUTION OF 
IONISED STRUCTURE 



Mesinger & Furlanetto (2007) introduced an approximate 



but efficient method for simulating the reionization process. 
This so-called semi-numerical method extends prior work by 
Bond & Myersl (119961) andlZahn et al.l (120071). The method 



generates an estimate of the ionization field based on a cat- 
alogue of sources assigned within the halo field by applying 
a filtering technique. Good agreement is found with numer- 
ical simulations, implying that semi-numerical models can 
be used to explore a large range of reionization scenarios. 
In this paper we apply a semi-numerical technique to find 
the ionization structure resulting from GALFORM galaxies 
within the Millennium-II dark matter simulation. 



3.1 HII regions 

We begin by binning galaxies from the GALFORM model 
into small regions of volume (or cells) . We assume the num- 
ber of photons produced by galaxies in the cell that enter 
the IGM and participate in reionization to be 



N. 



7, cell 



/, 



Jo 



Lye, cell 



(t) dt, 



(4) 



where / eS c is the escape fraction of photons produced by 
stars in a galaxy. The total Lyman continuum luminosity of 
the iVceii galaxies within the cell expressed as the emission 
rate of ionizing photons (i.e. units of photons/s) computed 
from GALFORM is 



iVi 



Lye, cell (t) — 2_^ ^L.yc,i(t), 



where 



Ni 



Lye 



i(t) = [° 



L U ,j(t) 

hv 



dv, 



(5) 



(6) 



L u ,i is the spectral energy distribution of galaxy i and ^thresh 
is the Lyman-limit frequency, /i^thresh = 13.6 eV. Note that 
the number of photons produced per baryon in long-lived 
stars and stellar remnants depends on the IMF and metal- 
icity (Z). The stellar population models used in GALFORM 
output give 4.54xl0 3 for Z=0.02 and 6.77xl0 3 for Z=0.004 
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Figure 1. The UV luminosity functions from the Bow06, Lagos and NOSN models with observed data points. The top-left panel shows 
the predicted UV luminosity functions from the models at z=10.074 with the observations estimated by |Bouwens et alTj pOll) (circle and 
1-cr upper limit arrows, 1600 A ) for z ~ 10. The top -right panel shows the z=7.883 predicti ons fr om the models including the observations 



for the 
predictions for z 



8 measured by Bouwens et al. 



( 2010) (triangles, 1700A) and |McLure et al.| (|2Q10|) (cir cles, 1500A). Th e bottom-left is the 
7.27 2 from mo dels with the observations for z ~ 7 came from |McLure et al.| ( |2010 > (circles, 1500A), Oesch et al. (2010) 
(triangles, 1600 A) and|Ouchi et al. (2010) (squares, 150 0A). The bottom-right is for z=6.197 predictions with the observatrions for z ~ 6 



measured by McLure et al. (2009) (circles, 1500A), and Bouwens et al. (2007| (triangles, 1350A) 



using the Kennicutt IMF. Note that we assume the total 
Lyman continuum luminosity in a cell at redshift z% to be 
constant until the next snapshot at redshift Zi+i, and calcu- 
late the number of photons produced in the cell between z% 
and z i+ i to be iV L yc,ceii(^zJ x (t z . +1 -t z .). 

We then calculate the ionization fraction within each 
cell according to 



Qcell = 



iV~ 



(7) 



_(l + F c )iVHI,cell. 

where F c denotes the mean number of recombinations per 
hydrogen atom up to reionization and iVHi, C eii is the number 
of neutral hydrogen atoms within a cell. The latter quantity 
is calculated as 

-7VHI,cell = ^Hl(^DM,cell + l)V^ e llj (8) 

where we assume that the overdensity of neutral hydrogen 



follows the dark matter (computed based on the Millennium- 
II simulation density field), nm is the mean comoving num- 
ber density of hydrogen atoms, and V^eii is the comoving 
volume of the cell. Self-reionization of a cell occurs when 
Qceii = 1- We divide the Millennium-II simulation box into 
256 3 cells, yielding cell side lengths of 0.3906/i _1 Mpc and 
comoving volumes of 0.0596/i _3 Mpc 3 . 

Theoretical prediction of the parameters F c and / eS c in 
Equations Q and (I7| is complicated, and their values are 
not known. The recombination parameter F c is related to 
the density of the IGM on small scales, while / esc depends 
on the details of the high redshift ISM. Previous work using 
GALFORM suggested the value (1 + F c )/f esc - 10 ( [Bern] 



son et al. 2001 Raicevic et al. 2011) to fit observational 



constraints on reionization. Here, we determine the value of 
(1 + F c )/ i 'esc that is required in order to give a particular 
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Figure 2. Maps of the 21 cm intensity in slices for a range of values of (xi) corresponding to different stages of reionization. We assume 
the Bow06 model. The units of the grey-scale are (28[(1 + z)/10]mK). The slices are 0.3906/i _1 Mpc deep. 



value of ionization fraction at each redshift. We explicitly 
note here our assumption that values of (l + F c )// esc do not 
depend on galaxy mass. In reality the escape fraction may be 
mass dependent, particularly in models with SNe feedback. 
Based on equation (J^TJ) , individual cells can have Q C eii > 
1. On the other hand, cells with Q ce ii < 1 may be ionized by 
photons produced in a neighbouring cell. In order to find the 
extent of ionized regions we therefore filter the Q ce ii field us- 
ing a sequence of real space top hat filters of radius R (with 
0.3906 < R < 100/i -1 Mpc), producing one smoothed ion- 
ization field Qr per radius. At each point in the simulation 
box we find the largest R for which the filtered ionization 
field is greater that unity (i.e. ionized with Qr > 1). All 
points within the radius R around this point are considered 
ionized. This procedure forms the position dependent ion- 
ization fraction ^ Q ^ 1, which describes the ionization 
structure of the IGM during reionization. The filtering fol- 
lows the method outlined in more detail in |Geil &; Wyithe] 
(2008). 



3.2 redshifted 21cm intensity 

Fluctuations in 21 cm intensity (or brightness temperature) 
from different regions of the IGM include contributions from 
a range of different physical properties, including density, 
velocity gradients, gas temperature, hydrogen spin temper- 
ature and ionization state (Furlanetto et al.||2006|). In this 



paper we restrict our attention to analyses that assume the 
spin temperature of hydrogen is coupled to the kinetic tem- 
perature of an IGM that has been heated well above the 



CMB temperature (i.e. T s ^> Tcmb). This condition should 
hold during the later stages of the reionization era (z < 9, 
|Santos et al.|2 007). In this regime there is a proportionality 
between the ionization fraction and 21 cm intensity, and the 
21 cm brightness temperature contrast may be written 



AT = 23i 



1 + 3 

10 



[1 - Q] (1 + <Weii) mK. (9) 



Here we have ignored the contribution to the amplitude from 
velocity gradients, and assumed as before that the hydrogen 
overdensity follows the dark matter (foM.ceii). 



4 STRUCTURE OF REIONIZATION 

In this section we present results for the possible structure 
of ionisation in the IGM. Figures [2] - [6j we show example 
ionization maps for our five different models. In each case 
we show examples for (xi) = 0.95, 0.75, 0.55, 0.36, 0.16 and 
0.056, illustrating the growth of HII regions during reioniza- 
tion. For our model these values correspond to redshifts of 
z - 6.197, 6.712, 7.272, 7.883 8.550 and 9.278 [selected for 
comparison with the work by Lidz et al7| (|2008[)]. Maps are 



shown for Bow06, NOSN and Lagos models in Figure [2] - [3] 
and for Bow06(no suppression) and NOSN (no suppression) 
models in Figure [5] -|6] In order to make this comparison we 
adjust the quantity (1 + F c )/ f esc in equations Q and ff\ so 
as to get same mass averaged ionization fraction (xi) for all 
models at each redshift. The values of (1 + F c )/ f esc required 
in order to give a particular values of ionization fraction at 
each redshift are shown in Table [2] Models presented in this 
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Figure 3. Maps of the 21 cm intensity in slices for a range of values of (xi) corresponding to different stages of reionization. We assume 
the NOSN model. The units of the grey-scale are (28[(1 + z)/10]mK). The slices are 0.3906fr -1 Mpc deep. 
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Figure 4. Maps of the 21 cm intensity in slices for a range of values of (xi) corresponding to different stages of reionization. We assume 
the Lagos model. The units of the grey-scale are (28[(1 + z)/10]mK). The slices are 0.3906/i -1 Mpc deep. 
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Figure 5. Ionization maps for a range of values of (xi) corresponding to different stages of reionization. We assume the Bow06(no 
suppression) model. The units of the grey-scale are (28[(1 + z)/10]mK). The slices are 0.3906/i -1 Mpc deep. 
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Figure 6. Ionization maps for a range of values of (xi) corresponding to different stages of reionization. We assume the NOSN (no 
suppression) model. The units of the grey-scale are (28[(1 + z)/10]mK). The slices are 0.3906/i -1 Mpc deep. 



© 2011 RAS, MNRAS 000,n]{T4 



Reionization in Galaxy Formation Models 9 



fiducial model 
Bow06 v cut = 30 km/s <Xj>=0.55 



fiducial - (NoSN, v cu{ = 30 km/s) 



fiducial - (NoSN, no supression) 



I 






(fiducial - Bow06, no suppression) 


40 


i ' i ■ i 


|0 


- 


-40 





(NoSN, v cut = 30 km/s) - (NoSN, v clA = km/s) 



Q. 



-40 




Mpc 



40 




-40 





fiducial - (Lagos, 


u cut = 30 km/s) 


40 


i ■ i 


^ ' ' 


|0 






-40 


i . i 


>> . 



-40 




Mpc 



40 



Figure 7. Example maps of the ionization structure produced by our modelling. In each case the slices shown are 100h -1 Mpc on a side, 
and 0.3906h -1 Mpc deep. The models were computed with (pa) = 0.55 (z = 7.272). The top-left panel shows map for the Bow06 model. 
To make the effect of SNe and photo-ionization feedback on the ionization map clear, the others show the subtraction of the Bow06 map 
from the variant model's maps except the bottom-middle panel. Here a positive value shows regions where the Bow06 model predicts 
HII region but variant models does not. Conversely the negative values show regions where the variant models predict ionization but the 
Bow06 model does not. The units of the grey-scale are (28[(1 + z)/10]mK). 



Table 2. The values of (1 + F c )/ f esc corresponding to the different models and redshifts shown in this paper. 



Redshift (z) 


9.278 


8.550 


7.883 


7.272 


6.712 


6.197 


{xi) 


0.056 


0.16 


0.36 


0.55 


0.75 


0.95 


Bow06 


4.85 


3.24 


2.72 


3.10 


3.83 


4.82 


Bow06(no suppression) 


4.856 


3.24 


2.71 


3.12 


3.85 


4.81 


Lagos 


3.86 


2.61 


2.17 


2.53 


3.11 


3.95 


NOSN 


417.98 


189.28 


106.70 


85.59 


74.83 


68.94 


NOSN(no suppression) 


267.78 


136.97 


85.03 


73.86 


69.62 


68.28 



work take values of (1 + F c )/f eS c that are less than 10 for 
models including SNe feedback in agreement with the work 



of Benson et al. ( 2001 ) and Raicevic et al. ( 2011 ), but greater 



than 50 for models without SNe feedback. 

The effects of SNe feedback strength (between the 
Bow06 model and the NOSN model) and of star formation 
law (between the Bow06 model and the Lagos model) on the 
ionization history can be seen explicitly by comparing Fig- 
ures [3] and [4] with Figure [2] The regulation of star formation 
and cooling of hot gas in small galaxies by the SNe feedback 
process leads to massive galaxies which are more biased to- 
wards dense regions, dominating the production of ionizing 
photons. As a result, the evolution of large HII regions in 
the Bow06 model (Figures [2| starts from the overdense envi- 



ronment and propagates to neighbouring overdense regions. 
Conversely, the production of ionization photons from the 
massive galaxies is less prominent in the NOSN model (Fig- 
ures p| than in the Bow06 model. As a result the HII re- 
gion evolution maps for the NOSN model show many more 
smaller HII regions than in the Bow06 case. The different 
star formation law prescriptions in the Lagos (Figures B| 
and Bow06 models also lead to differences in evolution of 
HII regions. However the variation is much smaller than is 
found from differences in the SNe feedback strength. 

The effect of photo-ionization feedback on the HII re- 
gions evolution can be seen by comparing Figures [5] and [6] 
with Figure |2] We find very little difference between the HII 
region evolution in the Bow06 model and the Bow06(no sup- 
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Table 3. The values of F^ hotons , the fraction of ionizing photons 
where the galaxy circular velocity less than 30km/s. 



Model 



Vhalo <30Km/s All 



Bow06(no suppression) 0.0028 1 

NOSN(no suppression) 0.39 1 



pression) model (Figure pi. Thus the effect of SNe feedback 
on the evolution of ionisation structure is much larger than 
that from photo-ionization feedback. However in the absence 
of SNe feedback, the effect of photo-ionization feedback ef- 
fect is significant (Figure [6]), and the NOSN(no suppression) 
model produces numerous, smaller HII regions that are rel- 
atively homogeneously distributed through the IGM. 

To highlight the differences between the maps produced 
by models with and without SNe and/or photo-ionization 
feedback, in Figure [7] we show differences between ionisation 
maps for the Bow06 and other models. The exception is 
the lower-middle panel, which shows the difference between 
maps for the NOSN and NOSN(no suppression) models. In 
all cases the mass averaged ionization fraction is (xi) = 0.55 
(z = 7.272). Positive values represent the area where the 
Bow06 model predicts HII regions but alternative models do 
not (in the lower-middle panel the positive values represent 
the area where the NOSN model predicts HII regions but the 
NOSN(no suppression) does not). This figure clearly shows 
the large effect of SNe feedback relative to photo-ionization 
feedback and star formation prescription. 



5 CONTRIBUTIONS TO THE IONISING 
PHOTON BUDGET 

Radiative feedback has been thought to play a significant 
role in self-regulating the reionization process by suppress- 
ing galaxy formation in reionized regions (e.g. Iliev et al. 
2007 ) . These studies were based on the assumption that the 



ionising luminosity to halo mass ratio does not depend on 
halo mass. However the presence of SNe feedback in galaxy 
formation models is known to modify the mass-to-light ratio 
of galaxies through regulation of star formation. Motivated 
by the unexpectedly small difference in ionising structure be- 
tween models with SNe feedback that do and do not include 
radiative feedback, in this section we calculate the effect of 
radiative feedback on stellar mass and ionising photon con- 
tribution. 

First, in Figure [§] we show the cumulative fraction of 
stellar mass as a function of Vhaio [km/s] at (xi) = 0.55 for 
the Bow06(no suppression) and the NOSN (no suppression) 
models. The contribution to total stellar mass from galaxies 
which have circular velocities Vhaio smaller than 30km/s (left 
dotted line) is almost zero for the Bow06(no suppression) 
model and 25 percent for the NOSN (no suppression) model. 
This shows that SNe feedback greatly lowers the potential 
contribution of low circular velocity galaxies, which are the 
ones affected by the photo-ionization feedback process, and 
explains why SNe feedback is the more dominant effect. 

Similarly, we also calculate the fraction Fn photons °f 
ionising photons produced by galaxies with Vhaio less than 



30km/s for the Bow06(no suppression) and the NOSN (no 
suppression) models. The values are listed in Table [3] and 
show the small fraction of photons produced by the low mass 
haloes which are subject to the radiative feedback process; 
see also Raicevic et al. (2011 ) for similar discussion. Thus we 



find that SNe feedback renders the effect of radiative feed- 
back on the reionization history negligible, indicating that 
reionization is not self regulating, as indicated in previous 
work (e.g. |Iliev et al.|2007| ). 

While the calculations presented above provide the 
quantitative estimate of the effect of SNe feedback on the 
photon budget, we can also provide a simple argument to 
show qualitatively why SNe feedback should be the domi- 
nant process governing the contribution of low mass galaxies 
to reionization (e.g Benson et al.|2006 Raicevic et al.|2011 ). 
Ignoring photo-ionization feedback (Vc U t=0), gas cooling is 
very efficient in halos with a virial temperature of T v i r ~ 10 4 - 
10 5 K. We assume that the resulting star formation timescale 
is shorter than the Hubble time, in which case all of the 
cooled gas will either form stars or be ejected by SNe feed- 
back (equation [2} , yielding 



(1-R + P)M* -M b , 



(10) 



where M* is the mass of stars formed (before recycling), 
< R < 1 is the recycled fraction, and Mb is the to- 
tal mass of baryons in a halo of mass Mhaio- For low 
mass galaxies (with V cu t ~ 30km/s) f3 ^> 1, and the frac- 
tion of baryons converted into stars is M*/Mb ~ 1//3 ~ 
10 -3 (Vhaio/30kms -1 ) 3 oc Mhaio- The total ionising con- 
tribution to reionization is proportional to the product 
of this fraction and the mass in dark matter halos (i.e. 
(M* /Mb) Mhaio oc M^aio). At low masses the halo mass func- 
tion (number density per unit mass) is dn/dMh a io oc M^ lo 
with 7^2. The mass in stars per logarithm of halo mass 
per unit volume in the Universe is therefore proportional to 
Mh a i x Mhaiodn/dMhaio oc M ha io. Thus, we find that very 
low mass galaxies should contribute little to reionization (see 
also |Wyithe fe Loeb|2012 ) . 



6 THE 21-CM POWER SPECTRUM 

The filtering procedure described in § [3] provides a 3- 
dimensional map of the ionization structure within the 
Millennium-II Simulation box, which provides a 3 dimen- 
sional 21cm intensity cube via equation Nn. From this cube 
we calculate the dimensionless 21-cm power spectrum 



A 2 (k) = k 3 /(27r 2 )P 2 i(k) 



(11) 



as a function of spatial frequency k, where P21 (k) is the 21cm 
power spectrum. When calculating the power spectrum, ve- 
locity gradients increase the amplitude of the spherically 
averaged redshift space power spectrum by a factor of 4/3 
relative to the real space power spectrum based on linear 
theory ( |Barkana fe Loeb||2005] ). |Mao et al.j fl2012| shows 
that this factor can be much higher than 4/3 over the in- 
termediate range k ~ 0.1 — 1/i/Mpc at the epoch where the 
IGM is 50% ionized. 

The results of Figures [2J6] indicate that the 21cm power 
spectrum will depend on the galaxy formation model as- 
sumed. This is shown in Figure [5] which displays power spec- 
tra for each of the semi-analytic models in Table^ From this 
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Figure 9. The predicted 21-cm dimensionless power spectra for the models discussed in this paper. Panels are shown for a range of 
values of (xi) corresponding to different stages of reionization as shown in Figures [2] -Im The units of the dimensionless power spectrum 
are (28[(1 + z)/10]mK) 2 . In the lower sub-panels we show the ratio of the the other models to the Bow06 model. 



figure we see that Bow06 and NOSN models show a large 
variation in 21-cm power spectrum predictions, with the am- 
plitude of 21cm power spectrum for the Bow06 model being 
higher than the NOSN model across all wave- numbers (with 
the exception of very late in the reionization process). Sec- 
ondly, the Bow06 and Lagos models show slightly different 
21-cm power spectrum predictions. This is because the mod- 
ified star formation law in the Lagos model relative to the 
Bow06 model leads to different predictions for the number of 
luminous galaxies (Figure [l]) and hence the clustering of the 
the ionising source population. Thirdly, the NOSN model 
has a larger amplitude for the 21cm power spectrum than 
does the NOSN (no suppression) model. This shows that the 
photo- ionization effect on the 21cm power spectrum can be 
seen in the no SNe feedback models. Conversely, we find neg- 
ligible difference between the Bow06 and Bow06(no suppres- 
sion) models. This very small difference means that photo- 
ionization feedback can only effect the reionization signature 
in the absence of SNe feedback. These findings represent the 
main results of the paper. 

Figure [5] also shows the evolution of predicted 21-cm 
dimensionless power spectra for the different models. The 
evolution of the power spectrum in Bow06, Bow06(no sup- 
pression) and Lagos models show the characteristic rise and 
fall described in detail by Lidz et al. (2008). The maxi- 



mum amplitude of the 21-cm power spectrum occurs at a 
scale of around k ~ 0.2/i _1 Mpc for an ionization fraction of 
(xi) rsj 0.75. In all models there is a trend for the wavenum- 
ber k at which the shoulder due HII regions occurs to de- 



crease (corresponding to increasing size of HII regions) with 
increasing ionisation fraction. The largest difference is seen 
between the Bow06 and NOSN models, and is most pro- 
nounced at large scales (i.e. small wave numbers). In this 
regime the NOSN power spectrum is lower than for Bow06. 
At smaller wavenumbers (i.e. large scales) the Bow06 model 
has higher amplitude than the NOSN model for all ioniza- 
tion fraction ranges. The difference in amplitude between 
the two models increases from (xi) ~ 0.056 to (xi) ~ 0.16, 
before decreasing later in the reionization history. We can- 
not distinguish the difference between the Bow06 and the 
Bow06(no suppression) model at any redshifts. There is a 
small difference between the Bow06 and Lagos models at 
all redshifts, although the magnitude of difference is smaller 
than between the Bow06 and NOSN models. There is also 
a difference between the NOSN and NOSN (no suppression) 
models. The bottom panels in Figure [9] show the ratio be- 
tween the models and the Bow06 model for each redshift. 



6.1 Observational implications 



Lidz et al. (2008) demonstrated that first generation low 
frequency arrays like the MWA[j should have sufficient sen- 
sitivity to measure the amplitude and slope of the 21 cm 
power spectrum. To quantify the effect of SNe feedback on 
the power spectrum we therefore compare the amplitude 
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Figure 8. The cumulative fraction of stellar mass as a function 
of Vhalo[km/s] at (xi) = 0.55 (corresponding to z = 7.272) for the 
NOSN(no suppression) and Bow06(no suppression) models. The 
vertical line indicates the value of V cu t = 30km/s. 
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Figure 10. Plots of the evolution in dimensionless 21cm power 
spectrum amplitude (lower panels) and slope (upper panels) as 
a function of ionization fraction (xi). Predictions are shown for 
models, Bow06 (triangles, black solid line), NOSN (pentagons, 
violet dotted line), Bow06(no suppression) (squares, blue long 
dashed line), Lagos (circles, red dashed line) and NOSN(no sup- 
pression) (octagons, green dot dashed line). Results are shown 
for two central wave numbers, k p = 0.4/i -1 Mpc (left) and 
0.2h~ 1 Mpc (right) , corresponding to the point on the power 
spectrum at which we evaluate the amplitude and gradient. 



and slope of predicted 21-cm power spectra for the Bow06, 
NOSN, Bow06(no suppression), NOSN(no suppression) and 
Lagos models. In Figure [To] we plot these values as a func- 
tion of the (xi) for central wave numbers of k p = 0.2/i _1 Mpc 
and 0.4/i _1 Mpc, corresponding to the range of wave numbers 
to be probed by the MWA. There are significant differences 
in the predicted quantities. For k p = 0.2/i _1 Mpc, the inclu- 
sion of SNe feedback results in fractional changes that are 
of order unity, particularly near the peak of reionization. 

Since the ionization fraction is not a direct observable, 
we plot the progression of a model in the observable plane 
of power spectrum amplitude vs slope. These are shown for 
the variant models in Figure [TT] again for the two values of 
central wavenumber k p , corresponding to the point on the 
power spectrum at which we evaluate the amplitude and 
gradient. The arrows show the direction from high to low 
(xi) (from 0.95 to 0.056). The tracks separate according to 
whether SNe feedback is included or not and difference of 
star formation prescription in the models (the Bow06 and 
Lagos models). To illustrate the potential for detectability 
of this difference we also include error bars at each point cor- 
responding to estimates for the MWA ( specifically an r~ 2 
distribution of 500 antennas) (Lidz et al. 2008| ) assuming 
1000 hours integration and 6MHz bandpasses for wavenum- 
ber /c p =0.4/i -1 Mpc. The figure demonstrates that mid-way 
through reionization (i.e. at the highest amplitude), the dif- 
ference between the tracks for models with and without SNe 
feedback could be detected by the MWA, indicating that 
the strength of SNe feedback during the epoch of reioniza- 
tion could be inferred directly from observations of the 21cm 
power spectrum. 



7 SUMMARY AND CONCLUSIONS 

Over the next decade we are likely to see the first measure- 
ments of the power spectrum of redshifted 21cm fluctua- 
tions from neutral hydrogen structure during the Epoch of 
Reionization. One goal of these experiments will be to learn 
about the properties of the galaxies that drove the reion- 
ization process. It is known that the ionisation structure of 
the IGM, and hence the observed 21cm power spectrum will 
be sensitive to the astrophysical properties of the reioniz- 
ing galaxies. With this in mind, Barkana (2009) has sug- 



gested that analytic models of the power spectrum could be 
used to determine the astrophysics of the reionizing galax- 
ies, provided they are tuned to provide a sufficiently precise 
description through comparison with numerical simulations. 
However, previous analyses of the structure of reionization 
and the predicted power spectrum have used very simple 
prescriptions to relate ionizing luminosity to the underlying 
dark matter distribution. In this paper we have made a first 
attempt to connect the details of the ionisation structure 
and 21cm power spectrum with realistic models for galaxy 
formation by combining the GALFORM galaxy formation 
model implemented within the Millennium-II dark matter 
simulation with a semi-numerical scheme to describe the re- 
sulting ionization structure. While not a true calculation 
of radiative transfer, semi-numerical models are known to 
reproduce the main features of the 21cm power spectrum 
where reionization is driven by UV ionizing sources. Our 
model includes a single value for the escape fraction of ion- 
ising photons from galaxies, independent of halo mass. In 
reality the escape fraction may depend on mass in a way 
that could be degenerate with the effect of SNe feedback. 
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Figure 11. Plots of the loci of points in the parameter space 
of 21-cm power spectrum amplitude and slope. Loci are shown 
for each of Bow06 (triangles, black solid line), NOSN (pentagons, 
violet dotted line), Bow06(no suppression) (squares, blue long 
dashed line), Lagos (circles, red dashed line) and NOSN(no sup- 
pression) (octagons, green dot dashed line) models. Results are 
shown for two central wave numbers, k p = 0.2h~ 1 Mpc (top) and 
0.4/i -1 Mpc (bottom), corresponding to the point on the power 
spectrum at which we evaluate the amplitude and gradient. The 
error bars for 0.4/i -1 Mpc at each point on the Bow06 model corre- 
spond to estimates for the MWA ( specifically an r~ 2 distribution 
of 500 antennas) ( |Lidz et al.|2008| > with 1000 hours of integration 
and 6MHz of bandpass. 



We find that the details of galaxy formation are re- 
flected in differences in the structure of reionization. As a 
result, each of the assumed star- formation law, radiative 
feedback and SNe feedback are found to affect 21cm power- 
spectrum predictions. Our main result is that the details of 
SNe feedback are most important in modifying HII region 
evolution, and hence the slope and amplitude of the 21 cm 
power spectrum. We find that photo-ionization feedback also 
affects HII region evolution but only in the absence of SNe 
feedback. Thus, unless SNe feedback is ineffective in high 
redshift galaxies, the reionization process is not self regulat- 
ing as has been argued previously (e.g. Iliev et al.|20 07). This 
finding is consistent with the work of Raicevic et al. (2011) 
who studied the photon budget in the context of the global 
evolution of reionization. We find that measurements of the 
amplitude and slope of the 21cm power spectrum would be 
sufficient to determine the level at which SNe feedback op- 
erated in high redshift galaxies. 

In this work we have concentrated on the effects of SNe 
and radiative feedback which are relevant to the galaxies 
thought to dominate reionization and are accessible to semi- 
analytic models implemented within the Millennium-II sim- 
ulation, and we have restricted ourselves to the assumption 
that the escape fraction is not mass dependent. This study 
illustrates the important role that semi-analytic models can 
play in realistic simulation of the connection between ion- 
ization structure and the properties of the galactic sources 
responsible for reionization. Our paper is the first step in a 



program to determine how redshifted 21cm observations can 
be used to probe astrophysics of reionization. 
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APPENDIX A: MODELLING SPATIALLY 
DEPENDENT REIONIZATION FEEDBACK IN 
GALFORM 

This appendix describes how spatially dependent reion- 
ization feedback is implemented in GALFORM for the 
calculations in this paper. The steps in the modelling are 
as follows: 

Step (1) We first run GALFORM at high redshift (from 
redshift ^20) to find the first resolved HII region using the 
scheme described in 



this step. 



3.1 



We assume (1 + F c )/ f esc — 1 for 



Step (2) For snapshots at redshifts where the first HII 
region is identified, we then find those galaxies which are 
inside HII regions and subject to radiative feedback through 
regulation of cooling processes (i.e. galaxies with V < V^ut). 

Step (3) A table listing these galaxies is generated for 
each snapshot. This table identifies which galaxies should 
have star formation regulated by radiative feedback during 
subsequent evolution up to the next snapshot redshift. 

Step (4) GALFORM is then run to the next redshift 
snapshot, including regulation of cooling processes in the 
galaxies identified in the table at step (3). 

Step (5) Steps (2) to (4) are then repeated for all snapshots 
down to redshift 6 where reionization is finished. 
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